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wavepackets -a long time phenomena- in the one dimensional quartic oscillator (a Kerr type Hamil- 
tonian). Two semiclassical theories are considered: time-dependent WKB and Van Vleck propaga- 
tion. We show that both approaches describe with impressive accuracy the autocorrelation function 
and wavefunction up to times longer than the revival time. Moreover, in the Van Vleck approach, 
we can show analytically that the range of agreement extends to arbitrary long times. 
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I. INTRODUCTION 

From the early times of quantum theory there has been 
a lot of interest in the possibility of describing quan- 
tum phenomena using approximate theories that take 
advantage of classical information [l|, Qj. Not only be- 
cause these semiclassical approaches provide in general 
a deeper insight into the system's behavior, due to the 
fact that classical quantities are more intuitive, but also 
in some cases the semiclassical computation of quantum 
quantities is made easier. 

The semiclassical propagation of wavefunctions started 
with the seminal work of Van Vleck 0], where a semi- 
classical propagator was first introduced. Gutzwiller Q 
showed that this propagator is the stationary-phase ap- 
proximation of the Feymnan path integral, which re- 
sults in a sum over classical paths. At the same time 
Gutzwiller provided some corrections to Van Vleck's for- 
mula which are essential for long times. 

The propagation of wavefunctions using the Van Vleck- 
Gutzwiller propagator has fundamental problems that 
manifest themselves more sharply in systems with a clas- 
sically chaotic dynamics. The facts that, for long times, 
the integrations to be performed are usually highly os- 
cillatory, and therefore not amenable to numerical com- 
putation 0, and that the number of orbits becomes 
unmanageably large, has raised doubts about the long- 
time semiclassical accuracy. These problems were con- 
trolled in Ref. [5[, where a semiclassical method was ap- 
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plied to compute autocorrelation functions of Gaussian 
wavepackets. This approach was shown to work both for 
chaotic [H and regular systems @, 0] in long-time regimes 
but it has not been applied to calculate wavefunctions. 

More recently, in Ref. Q it was shown that Gaussian 
wavepackets can also be propagated by using the stan- 
dard time dependent WKB theory (TDWKB) @, pro- 
vided it is complemented with short-time methods. In 
principle this scheme is suitable to compute wavefunc- 
tions in the long-time regime. Although initially con- 
ceived for chaotic systems, the scheme also works for 
nonlinear intcgrable dynamics (see Sec. IV below). 

Among the semiclassical time-domain methods one 
should also highlight the "initial value representation" 
(IVR) of the propagator, especially suited for numerical 
implementations (see Ref. |10| for a review). 

In spite of the great advances in semiclassical theo- 
ries in the last decades, there are still important open 
questions concerning the range of validity of semiclassi- 
cal approximations and their accuracy in the description 
of subtle interference quantum effects such as, for exam- 
ple, quantum revivals. 

The phenomenon of quantum revivals has been vastly 
studied in the literature (see Ref. [ll[ for a review) and 
observed in many experiments, from atomic and molec- 
ular to optical systems [Hj]. The Wigner function of an 
initially well-localized wavepacket spreads for short times 
in a classical way, then it enters a delocalized quantum 
regime, and eventually it recombines itself to recover its 
original form. Such a revival occurs at a time, T rev , 
which is long as compared with classical timescales, like 
oscillation periods. Perhaps more interestingly, in a wide 
class of circumstances, at times equal to a fraction of the 
revival time (pT rev /q) the wavepacket relocalizes into a 
number of smaller copies of the initial packet, giving rise 
to "fractional revivals" . When the initial wavepacket 
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can be associated to an essentially classical state, the 
"fractional revival" occurrence corresponds to a dynami- 
cal generation of Schrodinger's catlike states (a quantum 
superposition of macroscopically distinguishable states). 

One of the most interesting systems exhibiting both 
full and fractional revivals is the quartic oscillator, whose 
Hamiltonian is obtained by the squaring the harmonic 
oscillator Hamiltonian. In essence this is the Hamilto- 
nian of the Kerr model, describing a single mode of the 
quantized radiation field in a nonlinear medium, and ex- 
tensively studied in quantum optics. The formation of re- 
vivals and fractional revivals in this system was analyzed, 
for example, in Recently the quartic oscillator has 
experienced a renewed interest for its connections with 
quantum information processing in continuous variable 
(CV) systems. Indeed, Hamiltonians of the Kerr type 
are the simplest nonlinear ones acting on a single mode 
(i.e. one-mode quantum logic gate) needed to define uni- 
versal quantum computation within the subclass of uni- 
tary transformations generated by Hamiltonians that are 
polynomial functions of the CV operators [HI]. 

In addition to the already mentioned studies of revivals 
in the Coulomb problem @, 0|, we must also mention 
the papers by Wang and Heller [l5| and by Novaes 
The first authors considered the revival of a wavepacket 
in the Morse potential. Using a convenient numerical 
implementation of the Van Vleck-Gutzwiller propagator 
they succeeded in reproducing satisfactorily the first re- 
vival of the wavefunction. Novaes studied the semiclassi- 
cal propagation of a wavepacket in the quartic oscillator, 
starting from the semiclassical coherent-state represen- 
tation of the propagator. Even though he also obtained 
an excellent agreement, the increasing difficulty in deter- 
mining the required complex trajectories as time grows 
limited the application of the method to short times (a 
few classical periods) [l6[ . 

The present paper is devoted to show that "elemen- 
tary" semiclassical theories can be successfully applied 
to describe the revival phenomena in the quartic oscilla- 
tor. Our study focus both on the autocorrelation function 
and on the wavefunction. The two elementary semiclas- 
sical theories examined are two: Van Vleck propagation 
(Sec. Ill) and time-dependent WKB (Sec. IV). In the first 
case, calculations are analytical; in the second, numeri- 
cal. In both cases we find an excellent agreement between 
semiclassical theory and exact propagation even at very 
long times (e.g., multiples of the revival time). In the 
particular case of the Van Vleck autocorrelation function 
we prove analytically that it agrees with the exact one 
up to arbitrary long times, the error being semiclassically 
small and independent of time. As a byproduct of our 
study, we show that TDWKB also works efficiently in 
integrable nonlinear systems [8j]. 

Section II contains a description of the main aspects of 
the model. We present our main conclusions in Sec. V. 



II. THE MODEL 



Consider a one degree of freedom harmonic oscillator: 
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Throughout the paper we set m = 1 and u) = 1 (if de- 
sired, these constants can be recovered at any moment 
by dimensional considerations). We shall be interested 
in the dynamics of coherent states of this harmonic os- 
cillator, i.e., eigenstates of the annihilation operator [ij} 



(q + ip) , 



(2) 



with 



itl 



1 (a^ the creation operator). Our model to 
study the revivals and fractional revivals is the quartic 
Hamiltonian given by, 



H = 7 ?i 2 [ a)a + 



jh 2 



(3) 



where n = a^a is the number operator. This nonlinear 
Hamiltonian is of the Kerr type [Hj]. Recall that the 
effective Hamiltonian that describes the dynamics of a 
single light mode inside a high finesse optical cavity con- 
taining a Kerr medium is: 



H = uj\fi 



(4) 



(loi and uj' x are real frequencies), and we recover the op- 
tical context in our formalism considering Ti = 1. 

The quantum evolution of a coherent state, |«o), with 
the Hamiltonian (|3|) yields 



*Fn\ 



n=0 



(5) 

where \n) is a number state, i.e., n\n) — n\n). At multi- 
ples of the revival time 
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the dynamics reconstructs the initial coherent state. For 
times equal to (p/q)T2 (p/q an irreducible fraction) the 
evolved state consists of a superposition of q coherent 
states lying on a circle of radius |ao| (fractional revivals) 
[l3| . In Fig. [1] we show the Wigner function [lj| of the 
evolved state (O at several selected times. Schrodinger 
catlike states can be clearly seen at the fractional revival 
times. 

We shall focus on the autocorrelation function, 
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FIG. 1: (Color online) Snapshots of the evolution of an initially Gaussian wavepacket (Wigner functions). From left to right, 
times are t = T2/A, Is/3, T2/2. 38567, with T2 the revival time. The initial coherent state is centered at (qo,po) = (0, 14). 



in the semiclassical regime, defined by the condition 
|ao| 2 ^> 1, which is the appropriate semiclassical limit 
in the case of an optical context. It is important to note 
that this function shows the periodicity property 

C 1 (t + T 2 ) = e-^/ 4 C 1 (t) 7 (8) 

the period being the revival time ([6]). 

Note that the Weyl- Wigner representation 0, [2(| of 
the Hamiltonian © is: ~/{(q 2 +p 2 )/2} 2 - h 2 /4. Thus, the 
classical Hamiltonian corresponds to 

H(q,p)= 1 I 2 , (9) 

with the action variable 

I(q,p) = ^£ pdq=^(q 2 +p 2 ) , (10) 
a constant of motion. 



and the Maslov index /ifc, which (in the present case) 

coincides with the number of turning points (where q = 

0) encountered by the trajectory @. The amplitude is 
given by 

A k = ■ (13) 

with 

q" = q"(q',p',t). (14) 

Equation pip can be derived, for instance, by using time- 
dependent WKB theory to propagate a position eigen- 
state [1]. 

A. Autocorrelation function 



III. VAN VLECK APPROACH 

The semiclassical calculations in this section are based 
on the Van Vleck-Gutzwillcr approximation [3J, [J] to the 
propagator: 

p— i7r/4 

K(q", q', t) « V Ake iS k (Q',g',t)/h-i^/2 (n) 

k 

The sum runs over classical trajectories connecting q' to 
q" in time t. Each trajectory contributes with an ampli- 
tude Ak and a phase. The phase is made up from the 
Lagrangian action Sk, 

S k (q",q',t)= f (pq-H)dt, (12) 
Jo 



Consider the autocorrelation function 

/oc />oo 
/ dq"d q 'ro(q")Mq')K(q",q',t), 
OO J —OO 

(15) 

where ipoil) is the wavefunction of the initial state. By 
substituting the exact propagator K(q", q', t) by Van 
Vleck's we obtain a semiclassical correlation function. 

The simplicity of the system we are considering permits 
the analytical determination of the required trajectories; 
the properties of the trajectories can thus be calculated to 
the desired precision. The classical equations of motion 
are obtained from the Hamiltonian ©: 

q(t) = V2I sin[f3 + uj(I)t] , (16) 
p{t) = V2I cos[P + uj{I)t] , (17) 
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with 



c(/) = f = 2 7 I. 



(18) 



Both the action I -a constant of motion- and (3 are de- 
termined by the initial conditions. The trajectories con- 
tributing to the Van Vleck propagator are the solutions 
of the following boundary problem: 



q" = V2I sin[f] + uj(I)t], 
q' = V2Ism(f3), 



(19) 
(20) 



for a given time t. 

Among all trajectories that satisfy the equations 
above, only a small subset will be relevant to the au- 
tocorrelation function. The reason is that the initial 
wavepacket is localized on a phase-space region of radius 
0(\/~h), meaning that the important trajectories are only 
those having endpoints (q" ,p") and (q',p') in a region of 
radius O(Vh) around the center of the wavepacket. 

We shall consider, for convenience (but without loss 
of generality), an initial wavepacket corresponding to 
a coherent state |ao) located on the p-axis, i.e., ag — 
ipo/V2h, whose wavef unction is 



V'o(g) = (nti) 



-1/4, 



-q /2h e ip q/h 



(21) 



In this case, we can assume that both q" and q 1 are small 
[O(v^)] and resort to Taylor expansions around (go = 
0,po)- We can get a good idea of the structure of the set 
of trajectories that we need by analyzing the particular 
case q" = q' = 0. From Eqs. (|19I20[) we obtain a set of 
periodic trajectories, labeled by the winding number k, 
the corresponding actions being 



p 2 k kir 



The relevant values of k are those satisfying 
Ik ~ Iq 



Pi 



i.e., 



k 



2tt 



(22) 



(23) 



(24) 



For simplicity we ignore k — trajectories, which would 
require a special treatment. This means that our calcu- 
lation is not valid for times of the order or smaller than 
T\ , where T\ is the period of the classical motion of the 
centroid of the wavepacket: 



2lT _ 2lt 7T 
LOO 27/o 7/0 



(25) 



(This is not a problem as we are interested in long times.) 
The number of trajectories that contribute at t « kT\ is 



of the order of \k (see Sect. MI CI) . The properties of 
these periodic trajectories are easily calculated: 



2k. 
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It 



1 



(26) 
(27) 

(28) 



In the general case, when q" and/or q' are not strictly 
zero, even if the corresponding trajectories are not pe- 
riodic any more, they can be put into one-to-one corre- 
spondence with the periodic ones. Accordingly, the new 
actions and amplitudes can be calculated as Taylor ex- 
pansions in q" and q' around the periodic solutions. (The 
Maslov indices do not change, as they just count twice 
the number of turns.) 

We shall approximate the actions Sk to second order 
in (q" ',«/) and the amplitudes A% will be kept at zero-th 
order 0, 0). So, we need to calculate first and second 
derivatives of Sk with respect to q" and q' and evaluate 



(29) 
(30) 

, (31) 
(32) 



q" = 0,q' 


= 0,t). 


The results 


are 
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Thus, we arrive at 
S k 



n 2 k 2 
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The final steps in the calculation of the semiclassical cor- 
relation function are: (i) substitute the expressions above 
into Van Vleck propagator (fTT|) . (ii) insert the result- 
ing propagator into the definition of correlation function 
(fT5|) . (iii) calculate the Gaussian integrals. In this way 
we obtain: 



Ca(t) 
with 



-iir/A 



E' 

fe=i 
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(33) 



(34) 



and pk = W2kir/yt. 

Equation (|33| gives the semiclassical correlation func- 
tion obtained from Van Vleck's propagator and perturba- 
tion analysis around periodic orbits; it is a well behaved 
sum that can, in principle, be calculated for any t and 
compared with the exact result. The domain of validity 
is t ^ Ti. Now we proceed to numerical comparisons, 
postponing analytical considerations to Sect. IIII Cl 
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In Fig.[2l we display semiclassical and exact correlation 
functions in two time windows, one for short times [panel 

(a) ], the other for times around the second revival [panel 

(b) ]. Except for the initial correlation peak (t » 0), where 
the semiclassical approximation was expected to fail, the 
agreement is excellent. Remarkably, the peak that is lost 
at t w 0, it perfectly resurges at t w 2Ti (and at other 
revival times -graphics not shown). 




FIG. 2: Absolute value of the autocorrelation function vs. 
time, in units of the classical period T\. Line: exact; dots: 
semiclassical approximation. We set m,u,fr,y = 1, po = 
14. (a) Short times. Inset: blowup of the region where the 
first interference effects appear, (b) Times around the second 
revival (t « 2T 2 ). 



B. The Van Vleck wavefunction 

There is a straightforward test for the wavefunctions 
generated by the Van Vleck approach: If, instead of do- 
ing both the integrals defining the correlation function in 
(fTS")) . we calculate only the first one, over q' , we obtain a 
Van Vleck wavefunction ip(q"). First of all, this semiclas- 
sical wavefunction can be tested only in a small interval 
around q" — (because of the approximations we did 
in calculating the Van Vleck propagator). Moreover, this 
calculation omits the contributions arising from trajecto- 
ries that arrive at q" in time t with negative momentum. 
So, we expect the approximation to work well only at 
times where the (Wigner function of the) evolved state 
does not have negative momentum components around 



q" = 0. Such a situation arises, for instance, at re- 
vival times, when the initial wavepackets is exactly re- 
constructed. 

Another cases that can be described in terms of peri- 
odic orbits are fractional revivals like that at t — T2/3 
(see middle panel of Fig. [1]), when the wavepacket is 
a superposition of three coherent states with phases 
arg(a) = 7r/2, 2ir/3, 47r/3. In this case one has: 



W)) 



e~™' 12 co K) 



where 



3-27TI/3 



co 



(35) 



(36) 



and the ellipsis stand for the other two coherent states. 

Figure[3][(a) and (b)] shows that indeed the Van Vleck 
scheme which only uses periodic orbits reproduces the 
exact wavefunctions almost perfectly at times t = T 2 and 
i = T 2 /3. 

In order to calculate the wavefunction at an arbitrary 
time (but not too short), we must also take into account 
the contributions of the family of trajectories that, start- 
ing at q', arrive at q" in time t with negative momen- 
tum. As done before, we Taylor-expand actions and am- 
plitudes around the main family (q" = 0, q' = Q). These 
are "half-periodic" trajectories, i.e., their frequencies sat- 
isfy u>kt = 2n(k — 1/2), for k > 1. We skip the details 
and just show the final results: 
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, (37) 

(38) 
(39) 



Adding the contributions of the family above to that of 
the periodic trajectories we obtain a semiclassical wave- 
function valid for arbitrary times. Figure [3] (bottom 
panel) compares semiclassical and exact wavefunctions at 
a time when the state is completely delocalized in phase 
space (rightmost panel in Fig. [T]). Again, the semiclas- 
sical approximation performs remarkably well for small 
enough values of q. Not unexpectedly, for \q\ w 3 some 
small deviations start to show up, and keep growing with 
increasing 



C. Analytical comparisons 

The striking accuracy of the semiclassical approxima- 
tions verified in the previous sections, together with the 
relative simplicity of the expressions involved, suggests 
that it should be possible to give an analytical proof of 
the approximate quantum-scmiclassical equivalence. In 
the following we show analytically that the error commit- 
ted in the semiclassical autocorrelation function is really 
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FIG. 3: Real part of the wavefunction ip(q) vs. q at some 
selected times. Line: exact; dots: Van Vleck. Similar agree- 
ment is observed for the imaginary part. We used m = lo — 
% = 7 = 1, po = 14. (a) At the first revival time, (b) At the 
fractional revival time t = T2/3. (c) At t = T2/2. 38567, when 
the wavefunction is completely delocalized. 



very small and, what is especially important, indepen- 
dent of time. (A similar analysis could in principle be 
carried out for wavefunctions but will not be attempted 
here.) 

When comparing semiclassical (j3"3")l and exact ([7]) cor- 
relation functions we immediately see a fundamental dif- 
ference: time appears in a different way in both expres- 
sions. While the quantum expression is a Fourier series, 
i.e., a sum of plane waves (in t), the semiclassical correla- 
tion is a sum of wavepackets, or wavetrains. The tool that 
switches between wavetrains and plane waves is Poisson 
transformation fl9j . But before applying Poisson trans- 
formation to, say, the quantum correlation function let 



us introduce a small simplification. 



1. Quantum correlation function 

For large |ao| 2 one can approximate the Poisson dis- 
tribution by a Gaussian [T^ |: 



|ao| 2 ]£0 




-(n-y+1/2) 2 /2v 



(40) 



2 . With this approximation the exact 



where v = \a a 
autocorrelation can be written in terms of the Jacobi 
theta function [21| 



nin r-{-2inz 



(41) 



In fact, after the substitution Poisson-to-Gaussian in 
Eq. frj). we extend the lower limit of the summation to 
-co, obtaining 



d(t) 
where 



1 



-(.-1/2)72. e-^'/^Mn), (42) 
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2i 



1 - 
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2^ 



n = — ( -— - ijht 

m V Iv 



(43) 
(44) 



Expressing C\ (t) in terms of i9 3 brings in several ben- 
efits. First, formulas become more compact and can 
be calculated quickly and efficiently using standard soft- 
wares, e.g., Mathematica [22[, where $3 is a built-in func- 
tion. Second, Poisson transforming the theta function is 
equivalent to using the functional equation 



^(z| T ) = H r )-vv 2 /^ 3 (^|-i) 



(45) 



where (— it) -1 / 2 is to be interpreted by the convention 
I arg(— ir)\ < tt/2 [2l[. (This equation has already been 
used by Wang and Heller in their semiclassical study of 
the square well [H[.) 

Now we use the functional equation but restricting our- 
selves to times longer than the period T%, i.e., jKt 3> v. 
So, 



n 
1 

n 

z i 

{-ir t)- 1/2 



Zi 



IT ITT 

2 + 2jht + ' 
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1 -fht 2v 1 2 K 2 t 2 
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(46) 
(47) 
(48) 
(49) 
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In this way we arrive at an approximate expression for the 
quantum correlation function which has a semiclassical 
appearance: 



Ci(t) 



-i//2 



(50) 



Before comparing this expression with the semiclassical 
one we shall do some manipulations of the semiclassical 
formula. 



2. Van Vleck correlation function 

The semiclassical correlation function is written as a 
sum over periodic orbits. Each orbit is weighted by the 
exponential function: 



W(k) 



-(Pk-Po) 2 fa k H 



(51) 



Recall that this function depends on time through pk 
(|2"2")) . For long times (large k) we can set afc ~ 1. Thus 
W(k) becomes a real function, which can be well approx- 
imated by a Gaussian. 

We show in Fig. Q] some plots of the exact W(k) 
[Eq. (|5l"j) ] for typical values of the parameters together 
with the simplest Gaussian approximation obtained by 
setting k = fco, au = 1 in the denominator, and linear 
expanding pk around fc = fco> i-e., 



W(k) 



3 -7r(fc-fc ) /2fe 7fit 



(52) 



It is verified that the Gaussian approximation is already 




FIG. 4: Orbit weight function for times t/Ti = 3,5,7 (from 
left to right, respectively). Lines: Gaussian approximation, 
with k a continuous variable. Dots: absolute value of exact 
weight. Parameters are m,u},"/,h — 1, po = 14. 

very good for the small values of k considered. 

The remarkable fact is that if we substitute the approx- 
imate weight function ([52")) into the semiclassical correla- 
tion function (|33[) we obtain the approximate quantum 



correlation function ([50)1 . This concludes the analyti- 
cal comparison of quantum and semiclassical correlation 
functions: We have shown that for long times the differ- 
ence between both is semiclassically negligible. 



IV. TDWKB APPROACH 

It was recently shown that when the underlying clas- 
sical dynamics of a system is chaotic the evolution of a 
Gaussian wavepacket can be described by the standard 
time-dependent WKB (TDWKB) method [|. The rea- 
son is that, in semiclassical regimes, after some short time 
the Gaussian wavepacket stretches over a classical length 
~ h°. Thus, it becomes a primitive WKB state: 

M<i) = Mq)eMiSo(q)/h} , (53) 

supported by the Lagrangian manifold defined by 

P = dS /dq , (54) 

and where the amplitude Ao(q) and the phase So(q) are 
smooth function on the quantum scale. From then on the 
wavefunction evolves according to the TDWKB recipe: 



Mi) 



4%)exp^ 



(") 



(qyti-iHvir/2], (55) 



where v labels the different branches of the Lagrangian 
manifold obtained by evolving classically the initial man- 
ifold (PI) . 

In chaotic systems the applicability of TDWKB to 
wavepackets is guaranteed by the exponentially fast 
stretching of phase space [|[ . We show here that the same 
scheme can be applied in the case of integrable nonlinear 
systems, where the stretching is linear in time. 

For the present analysis we found more convenient to 
make a slight modification of the Hamiltonian of previous 
sections: 



H = ~fh 2 (n — uqY 



(56) 



The corresponding classical Hamiltonian is 

H(q,p)= 1 {I-I ) 2 . (57) 

This is equivalent to working with the Hamiltonian ([3]) 
but in the interaction representation, with a free evolu- 
tion given by the harmonic oscillator Hq = hui'n + C 
[with lu' — 2jh(no + 1/2) and C an appropriate con- 
stant)]. Thus, we eliminate the rotation dynamics of 
the wavepacket, while preserving the nonlinear squeez- 
ing. This choice simplifies the determination of the initial 
WKB manifold (|54[) . which now remains almost station- 
ary; otherwise it would rotate, forcing us to change repre- 
sentation from time to time to avoid caustics. Switching 
to the interaction representation does not affect the semi- 
classical accuracy of the calculation, because the trans- 
formation generated by Hq is semiclassically exact [23| ■ 
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The centroid (qo,Po) of the initial wavepacket will be 
chosen in such a way that I = h(n + 1/2) = (q +p )/2. 
The numerical implementation of the TDWKB recipe Q 
requires the determination of the initial manifold through 
Eq. ([54]) . The initial action So(q) is extracted from the 
phase of the exact wavefunction propagated up to some 
short time ti. The only condition this time must satisfy 
is that the exact wavefunction must be described to good 
accuracy by a primitive WKB state, i.e., Eq. (|53|) with 
A(q) and S(q) smooth on the quantum scale. 

Our choice of the initial manifold is showed in Fig^a). 
Using the classical equations of motion we evolve this ini- 
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FIG. 5: (a) (Color online) Wigner function of the evolved ini- 
tial coherent state, centered at the point (qo,Po) = (0, 14), for 
a time t = T2/32O (T2 the revival time). The superimposed 
black line is the initial Lagrangian manifold, (b) Wigner 
function of the same initial coherent state at t = T2/I6. The 
spiraling black line is the classical evolved manifold that sup- 
ports the WKB state. Parameters are m, 7, H = 1. See text 
for explanation of the meaning of crosses and circles. 

tial manifold up to the desired final time [see Fig. 0(b)]. 
In order to calculate the WKB wavefunction ipt(q) we 
have to determine the classical trajectories correspond- 
ing to each term in the sum in Eq. (|55[) . Those are 
the trajectories that at times U have initial conditions 
{q^f 1 > ) on the initial manifold, and at time t reach the 



final manifold at the points (qi ,Pj ), where q^ 1 1 = q 



If) _ 



for all values of j. For example, in Fig. ((5|) we show with 
crosses in panel (a) the initial conditions of the trajec- 
tories that end at points (q. 



(/) _ 



QiPj ) with positive 



momentum [panel (b)], and with circles the trajectories 
that end with a negative momentum. The classical ac- 
tion for these trajectories can be calculated analytically 
for the Hamiltonian in Eq. (157 



S. 



(3) _ 



So + J ( P q - H) dt = So + hlpWJp 

MlAlj - Hj)At , (58) 



where = 2j(Ij — Io), Hj = — Iq) 2 and At = t — 

U. The Maslov index fj, in Eq. ([55)l equals +/— the num- 
ber of turning points along a clockwise/counterclockwise 
trajectory @. The amplitudes in Eq. ([55|) are calcu- 
lated from the continuity equation 



U)rJf) 



A? (q 



= q)=A Q U i) ) 



dqi 
dq f 



1/2 



(59) 



where A (q^) is the amplitude of the primitive WKB 
wavefunction supported by the initial manifold. The 
factor \dqi/dqf\ is determined by evolving numerically 
a nearby trajectory. 

Note that the relatively hard part of this numerical 
method is the calculation of the points (qj ,Pj) and 

{q^ = q,Pj ) in the initial and the final Lagrangian 
manifolds, respectively. This can be done in a system- 
atic way defining a parameter s running along the initial 
manifold and calculating the points {q^\p^P) as inter- 
sections of the trajectories, with parameter s and fixed 
elapsed time t, with the phase-space vertical line q = , 
in a way resembling a Poincare section map. 

In Fig. [6] we compare the exact wavefunction of an ini- 
tially coherent state with the TDWKB wavefunction for 
t = 4X2 m panel (a) and for a generic time in panel (b) 
[when the Wigner function is nonlocalizcd, similar to the 
wavefunction in the rightmost panel of Fig. {TJ]. In panel 
(a) we see that the agreement is very good for all values 
of q, with small errors for \q\ > 1.8. For multiples of the 
revival time, t = mX^ the number of classical trajecto- 
ries needed to build up the TDWKB wavefunction grows 
like ex 90 m. Nevertheless, the numerical errors seem to 
be almost constant. Like in the Van Vleck approach of 
Sec. Ill, at multiples of the revival time only classical 
trajectories with positive final momenta contribute. 

In Fig. Elb) we show the TDWKB approximation in 
the generic case, where trajectories with final negative 
momenta also contribute. In this example, the exact 
wavefunction spreads also over the regions — 17 < q < 
— 11 and 11 < q < 17, where the final manifold crosses 
several times the axis p = 0. These points are caus- 
tics of the TDWKB approximation, where the amplitude 
in Eq. (|59p diverges. Thus, each time the contributing 
trajectories has a null final momentum, the TDWKB ap- 
proximation breaks down. This is clearly seen in Fig [S] 
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FIG. 6: Real part of the evolved wavefunction of an initially 
coherent state centered at the point (qo,po) = (0, 14) evolved 
with the Hamiltonian of Eq. (|56[) . The full line is the exact 
wavefunction and the circles correspond to the TDWKB ap- 
proximation, (a) For t = 4T2 (T2 is the revival time); (b) 
For T2/2. 38567. Similar agreement is observed for q < and 
for the imaginary parts. Parameters are m, 7, h = 1. 

(b). The proper treatment of the WKB function in 
this region requires more sophisticated approximations. 
For our present purposes, it is enough to verify that for 
|<7| < 10, where there are no caustic points, the agreement 
is excellent. 



V. CONCLUDING REMARKS 

Quantum revival is a subtle phenomenon where in- 
terference plays a crucial role. We studied revivals in 
the quartic oscillator from a semiclassical perspective. 
Among various semiclassical theories existing in the liter- 
ature we chose two of most basic and popular: Van Vleck 
propagation and time-dependent WKB. In both cases the 
results were impressive: quantum dynamics -in particu- 
lar, revivals- can be semiclassically described with great 
accuracy. 



In the TDWKB approach, we computed the wavefunc- 
tion numerically for times beyond the first revivals. Ex- 
cluding an initial stage (where standard WKB fails), the 
classical skeleton of the wavefunction appeared to be a 
spiraling manifold. Thus, we have exhibited another suc- 
cessful test of the TDWKB scheme for the propagation 
of wavepackets Q -in this occasion for the long-time dy- 
namics of an integrable system. However, at present, we 
cannot assess analytically how far the agreement will ex- 
tend. In order to do this one should integrate the present 
scheme with a theory capable of describing the short time 
dynamics. The natural tool is complex TDWKB which 
uses manifolds in the complexified phase space [24j . Then 
one should prove that, provided the dynamics is stretch- 
ing, the complex manifold describing a coherent state 
eventually decays into a real manifold [like that shown in 
Fig. 5(a)]. 

Like in some previous studies @, 0], we obtained an 
expression for the autocorrelation function as a sum over 
classical trajectories, all the ingredients being given in 
closed form. Furthermore, we showed analytically that 
the Van Vleck correlation function and the exact one are 
essentially equal. The key step was to use the Poisson 
transformation, which reshapes a semiclassical correla- 
tion function into a quantum-looking one, or vice versa. 
Similar analyses may be possible both for the Coulomb 
potential @, 0] and the Morse oscillator [l5[ (the latter 
is simpler, in principle, because its spectrum is quadratic 
in the quantum number [251]). 

The system under study, the quartic oscillator, is spe- 
cial, even among integrable systems, in that its Hamil- 
tonian is a (quadratic) function of the action variable. 
Thus, for instance, stationary WKB theory gives the ex- 
act energy levels for this system. We have provided an- 
alytical and numerical evidence showing that the semi- 
classical time- dependent schemes considered in this paper 
are also "exact" . It remains to ascertain if these consid- 
erations extend to more general H{I) Hamiltonians, e.g., 
of the polinomial type. 

Concerning the work by Novaes , who studied the 
quartic oscillator by using the semiclassical coherent- 
state representation of the propagator, it is likely, in the 
light of our results, that it should be possible to identify 
the relevant subset of complex trajectories which con- 
tribute to the autocorrelation function at long times. 
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